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Abstract 

Exact results of the finite-size behavior of the susceptibility in three- 
dimensional mean spherical model films under Dirichlet-Dirichlet, Dirichlet- 
Neumann and Neumann-Neumann boundary conditions are presented. 
The corresponding scaling functions are explicitly derived and their asymp- 
totics close to, above and below the bulk critical temperature T c are ob- 
tained. The results can be incorporated in the framework of the finite-size 
scaling theory where the exponent A characterizing the shift of the finite- 
size critical temperature with respect to T c is smaller than 1/v, with v 
being the critical exponent of the bulk correlation length. 

PACS: 64.60.-i, 64.60.Fr, 75.40.-s 



1 Introduction 

The basic ideas of the phcnomcnological finite-size scaling theory at criticality 
have been suggested by Fisher pQ, and Fisher and Barber 2 (for more recent 
reviews consult [3], @|, and According to the phenomenological theory, 

rounding and shifting of the anomalies in the thermodynamic functions set in 
when the bulk correlation length ^ becomes comparable to the characteristic 
linear size L of the system. More specifically, it is predicted that finite-size 
effects are controlled by the ratio £/£oo- Here we recall some fundamental 
notions and facts of that theory. 

*e-mail: daniel@imbm.bas.bg 
^e-mail: brankov@imbm.bas.bg 
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Let us start with a system having the geometry L d ~ d x oo d , with d! > di, 
where di is the lower critical dimension of the corresponding class of bulk sys- 
tems. Then, under boundary conditions r imposed across the finite dimensions 

of the system, the finite-size system exhibits a phase transition at a tempera- 
te) • • • (t) 
ture T = T c L , and the corresponding infinite system at T = T^do = T c . The 

so-called fractional shift, characterizing the shift of the critical temperature of 

the finite-size system is defined as 

e¥=(T c -T$)/T c cb^L-\ (1.1) 

where the expected asymptotic behavior for L 3> 1 is given by the shift exponent 
A. 

Let us now consider the susceptibility \ (per spin) which in the bulk (infinite) 
system has a critical-point divergence of the type 

Xoo {T)~Ar\ i^0+, (1.2) 

where 7 = 7(d) is the c?-dimensional critical exponent and t = (T — T c )/T c . 
On approaching the finite-size critical temperature from above at fixed L one 
should have 

x { Z\T)~A¥i-\ i^0+, (1.3) 

where i= (t — /T c = e ( [ ] +t and 7 = j(d') (in general, 7 ^ 7). Let t[ t) l 

denote the temperature at which the considered finite-size property \^ (T) first 
shows significant (of the relative order of unity) deviation from its bulk limit 
Xoo(T). Then one defines the fractional rounding 6^ = [t[ t) l - T c ) /T c ~ 

cS T ^L~ e , L > 1. The "rounding" measures the region of crossover from bulk 
d-dimensional to d'-dimensional critical behavior. 

The basic assertions of the "orthodox" phenomenological finite-size scaling 
are that (i) the only relevant variable on which the properties of the finite- 
size system depend in the neighborhood of T c is L/£oo(T) ~ Lt u , and (it) the 
rounding occurs when £ao{T) — L. 

It is easy to see that assumption (ii) leads directly to the conclusion that 
6 = l/i/ and, from (i), it immediately follows that 

X i T) (r)^L^lM(L/eoo(T)), (1.4) 

or, equivalently, 

X { l ) {T)~L-'l v X^\tL 1 / v ). (1.5) 

Here X^ T \x) is the universal finite-size scaling function describing the critical 
behavior of %, where, in order to reproduce the behavior described by Eqs. (|1.2f> 
and (|1.3fl one should have: 

X (t) (x) ~ Xoox -7 , as x-^00, (1.6) 
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and 

X {T \x)~X^ ] x-^ as x ->0+. (1.7) 

It has been considered, see 511313], that a more general formulation of the 
finite-size scaling hypothesis is given by the equation 

Xl\t) ~ L^X^ftL 1 ^). (1.8) 

Apart from the allowed shift of T c from T c> £, Eqs. l|1.5f) and (|1.8II are equivalent 
if £oo diverges algebraically with exponent v > 1/A. 

We emphasize that the use of the shifted temperature variable i in the above 

(t) 

finite-size scaling hypotheses allows for any L-dependence of the shift e L , i.e. 
the shift exponent A remains arbitrary. The assertion that the only criterion 
determining the finite-size scaling effects in the critical region is £,oo{T) ~ L 
leads to the equalities A = 9 = \ jv. This result follows from the renormalization 
group derivation of finite-size scaling 0,0, see also 0. Except in some special 
cases (ideal Bose gas and spherical model with a film geometry and Dirichlet- 
Dirichlet [S], [S], E3 or Neumann-Neumann boundaries ^U], when one 

has a logarithmic shift of the type ±\nL/L for d = 3 and A = 1 in all other 
dimensions d > 2), this relation seems to be quite generally valid. We stress, 
nevertheless, that the relationship A = \ jv is not 2 a necessary condition for 
the finite-size scaling to hold in general. 



Phenomenological finite-size scaling for systems with large 
critical shift 

Let us now consider in a bit more details what will be the consequences if the 
shift is asymptotically large, i.e. when \ jv > A. From Eq. i|1.8fl °ne then has 

Xl\t) ~ Ui' v X^{tL x l v + b {T ^L l / u - x ). (1.9) 

Obviously, in order to make explicit statements, we have to consider the two 
possibilities i) M r ) > and ii) < separately. 

Case i) : > 

From Eq. i|1.6|) it immediately follows that at T c 

L 7 \ (1.10) 

i.e., the divergence of the susceptibility at T c with respect to L will be reduced 
in comparison with the " standard" behavior 

xi T) (T c )~xM(0)L^ (1.11) 

predicted by Eq. I|1.5fl . 

An example of a model with large positive shift of the critical temperature 
is the spherical model under Neumann-Neumann boundary conditions |ll)j (see 



x„ 



b (r) 



3 



below). For d — 3 the shift in the dimensionless critical coupling is equal to 
lnX/(47rL). For such a shift one immediately obtains that 



x ( 1\t c ) 



4ttL 
InL 



(1.12) 



Case ii) : &( T ) < 

Then, at T c one has tL 1/u = 6 (r) J L 1/l/_A -» -oo when L — * oo. Obviously, in 
order to give a general answer what will be the behavior of the susceptibility in 
this case, one needs to know the asymptotics of the scaling function X^ T \x) for 
x — > — oo. 

Generally speaking, when x = tL x l v — > — oo takes place, the behavior of the 
zero-field susceptibility in a 0(n) model depends on the fact whether n = 1, 
or n > 1, and on the geometry of the system. Then, under periodic boundary 
conditions, summarizing the results of EH]j |14| . [T5] for the Ising type models 
and that ones of 0, [E|, Ej, [El, [IHj, [EH for O(n), n > 2 models, one has 



xi p) (T<T c )~ { 



L d , 

Lexp[L d 
r2(d-l) 

) 

L 2 exp[cL d 



l a{T)lk B T), 
2 T{T)/k B T], 



d! = 
rf' = 1, 
tf = 1, 
d' = 2, 
if rf' > 1 and 
or d' > 2 and 



1 

71 > 2 
n > 2 
n = 1, 
n > 2, 



(1.13) 

where ct(T) is the interfacial (or surface) tension in the Ising model, T(T) is its 
analog - the helicity modulus - for an 0(n), n > 2 system, and c is a constant. 
For other than periodic boundary conditions the situation is not so clear and 
a detailed information for the behavior of the zero-field finite-size susceptibility 
is still to be established, but one can hope that the leading-order behavior will 
remain unchanged (eventually, the power of L in front of the exponential terms 
may change when d! = di). 

We shall demonstrate the consequences of the shift on the finite-size behavior 
of the susceptibility for a fully finite and a film geometry only. The interested 
reader can easily complete the list of all possible geometries. 



The fully finite geometry 

In order to reproduce the size dependence of the susceptibility at T < T c , one 
should have X^(x) ~ X ( ff ] x dl, -~< . Then at T c one obtains 



x ( 1\t c ) 



X. 



j^d-\v(d-2+ri) 



(1.14) 



For the spherical model under Dirichlet boundary conditions, taking into ac- 
count that A = 1, f] = and v = l/(d — 2) (for 2 < d < 4) one derives 

X { ?{T c )~AL d -\ (1.15) 
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where A is a constant. This is exactly the result recently reported in |12| . 
The film geometry 

Similarly to the above, one can write down again the general expressions, but 
since they are rather cumbersome, we shall consider only the case of a film 
geometry on the example of the three-dimensional spherical model. Taking into 
account that then /3Y(T) = K — K c j5U], where K is the dimensionless 
coupling in the system, and that = — lnL/(47ri), c = Air [8], [Sj, \12\. 
we obtain 

X { l\t c ) ~ L 2 exp(clnL/47r) ~ L\ (1.16) 

which is exactly the result recently reported in ^21 ( see a ls° below) . Note that 
now x^i) diverges at T c faster than in the case of the "standard" finite-size 
scaling prediction given by Eq. lp~TTl (we recall that 7 = 2/(d - 2) for the 
spherical model, i.e. 7/V = 2). It is clear that when d! = 2, but d > 3, one will 
have Xl C^c) ~ L 2 exp(cL d ~ 3 ), where c is a constant, and we have taken into 
account that A = 1 for all d > 2. 

Before finishing this overview of the general form of the phenomenological 
finite-size scaling theory and its consequences, let us just remind the reader 
what the finite-size behavior in the case of geometry L d ~ d x oo d with d! < di is. 
Then, in the finite-size system there is no phase transition of its own, therefore, 
7 = 0. The hypotheses stated above still hold, provided we set 7 = and 

(t) (t) 

replace T c )£ by appropriately defined pseudocritical temperature T^ L . The 
latter can be defined, for example, as the temperature at which the susceptibility 
reaches its maximum. Then, in the case of an algebraic bulk singularity of the 
type D~2J), one has, e.g., Xl\t^ l ) ~ X^Uil", L > 1. This asymptotic 
behavior is exploited in one of the basic methods for evaluation of bulk critical 
exponents from finite-size data. Here we would like to stress that T^ L depends 
on which physical quantity has been used for its definition, since the maxima of 
the susceptibility and, say, the specific heat take place at, generally speaking, 
different temperatures. 

In the next sections we will present exact results for the finite-size behavior 
of the zero-field susceptibility of the mean spherical model with a film geometry 
when Dirichlet-Dirichlet, Neumann-Dirichlet and Neumann-Neumann boundary 
conditions are applied at the surfaces of the film. The corresponding scaling 
functions will be derived and their asymptotic behavior will be analyzed. 
We will demonstrate the important role of the shift when it is asymptotically 
larger than L~ x l v . 

The structure of the article is as follows. In Section [21 we briefly present 
the model and the basic facts needed for the analytical treatment of the free 
energy and the susceptibility for a fully finite system. Section [3] contains the 
corresponding modifications of these expressions in the case of a film geometry. 
The finite-size critical behavior of the susceptibility is analyzed in Section ^ 
The article closes with a discussion of the results obtained and their eventual 



5 



generalization. Some technical details needed for evaluation of different sums 
are described in the Appendix. 

2 The model 

We consider a d-dimensional mean spherical model with nearest-neighbor fer- 
romagnetic interactions on a simple cubic lattice. At each lattice site r — 
(n, r 2 , • • • , r d ) € Z d there is a random (spin) variable a(f) 6 R and a A = 
{cr(r),r <E A} is the configuration in a finite region A C Z d , containing |A| sites. 
The boundary conditions (to be denoted by the superscript r) define the inter- 
action of the spins in the region A with a specified configuration {cr(r), r G A c } 
in the complement A c = Z d \ A. In the remainder we take A to be the paral- 
lelepiped A = C\ x £ 2 x • • • x with Ci = {1, . . . , L{\, and explicitly study 
the case of film geometry which results in the limit L 2 , • • • , L d — * oo at finite 
values of Li = L. In the finite r% direction it suffices to specify the values of 
er(0,r 2 ,- • ■ ,r d ) and a{L + 1, r 2 , • • • , r d ) for all (r 2 ,---,r d ) G C 2 X • • • x £ d . 

The finite-size scaling behavior of the mean spherical model has been studied 
so far under periodic, antiperiodic, Dirichlet, 0, [H], [2D> E2> Neumann [TO] . 
|2*H] and Neumann-Dirichlct IJj] boundary conditions (for a review see, e.g., [H])- 
For lattice systems under Dirichlet boundary conditions we mean that 

<J(0, r a , • • • , r d ) = a(L + 1, r 2 , ■ ■ ■ , r d ) = 0, (2.1) 

under Neumann boundary conditions that 

0-(O, r 2) • • • , r d ) = <t(1, r 2 , ■ ■ ■ , r d ), a(L+l, r 2 , ■ • • , r d ) = a(L, r 2 , • • • , r d ), (2.2) 

and, under Neumann-Dirichlet boundary conditions that 

cr(0, r 2 ,---,r d ) = er(l,r 2 , • • • ,r d ), a(L + 1, r 2 , • • • , r d ) = 0. (2.3) 

Obviously, the terminology Dirichlet and Neumann boundary conditions is jus- 
tified by analogy with the continuum limit. The case of free surfaces in a system 
of film geometry (in the limit L 2 , — , L d — > oo), considered in the literature [S], 
|H] i [2D; [22]; corresponds to Dirichlet boundary conditions (r = D). When a 
fully finite system is envisaged, we always assume periodic boundary conditions 
with respect to the coordinates r 2 ,- ■ ■ ,r d , i.e., for all r G A and all integers 
m 2 , • • • ,m d , we set a(ri,r 2 + m 2 i 2 , • • • ,r d + m d L d ) = a{r x ,r 2 , ■ ■ ■ ,r d ). 

For brevity of notation, we consider the configuration space fl A = R' A ' as an 
Euclidean vector space in which each configuration is represented by a column- 
vector <t a with components labelled according to the lexicographic order of the 
set {(ri,r 2 , • • • ,r d ) G A}. Let <r+ be the corresponding transposed row-vector 
and let the dot (•) denote matrix multiplication. Then, for given boundary 
conditions r = (t%,t 2 , - ■ ■ ,T d ), specified for each pair of opposite faces of A 
by some r, = p (periodic), D (Dirichlet), N (Neumann), or ND (Neumann- 
Dirichlet), and given external magnetic field configuration h A = {h(r),f G A}, 
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with h(r) 6 R, the Hamiltonian of the model takes the form 



f3H ( A } (a A \K, h A ; s) = -\Ka\ ■ ■ a A + s o{ ■ a A - h{ ■ a A . (2.4) 

Here p = l/fc^T is the inverse temperature; K = [3 J is the dimensionless 
coupling constant; s is the spherical field which is to be determined from the 
mean spherical constraint, see equation l|2.19|l below; the |A| x |A| interaction 

(t) 

matrix Q A can be written as 



(A[ Tl) +2E 1 ) x ... x (A d T "'+2E d ), 



(Td) 



(2.5) 



where x denotes the outer product of the corresponding matrices, A I - T ''* is the 
Li x Li discrete Laplacian under boundary conditions r^, and Ei is the Li x Li 
unit matrix. 

As it is well known, the complete set of orthonormal eigenfunctions, {uy'(r, k), 
k = 1, . . . , L}, of the one-dimensional discrete Laplacian for periodic, Dirichlet, 
Neumann and Neumann-Dirichlet boundary conditions is given by 



,(p) 



(r, k) = IT 1/2 exp [-i rip*?) (k) 



,(N) 



and 



(r,fc) = 



AND) 



i[ D \r,k) = [2/(L + l)]V 2 sm \np[ D \k) 

L- 1 ' 2 for k = 1 

(2/L) 1 /2 C0S [( r _ i)^W(fc)] for k = 2, • 



L 



(r,fe) = 2{2L+ l)" 1/2 cos 



r- 5 



(ND) 



(A;) 



where 



yf)(fc) = 7r(fc ~ 1) 



7r(2/c - 1) 



L 1 TL w 2L+1 
The corresponding eigenvalues are 

\[ T \k) = -2 + 2cos^ r) (fc), fc = 1, . . . , L. 

The eigenfunctions of the interaction matrix Ij2.5|l have the form 



u[ T) (r, k) = u^XnM) u { ^\r 2 , k 2 ) 
and the corresponding eigenvalues are 

d 



A T d): 



r<i,k d ), 



{T) (k) = 2 ^ cos ^ (K), fee A. 



(2.6) 
(2.7) 
(2.8) 

(2.9) 

(2-10) 
(2.11) 

(2.12) 



k e A, (2.13) 



(2.14) 
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In order to ensure positivity of all the eigenvalues — ^Kfjff'Qc) + s, k € A, 
of the quadratic form in f3'H A T \a A \K, h A ;s), see equation 1)2. 4[l . it is convenient 
to introduce a shifted and rescaled spherical field > by setting 

8 = 8(4):= ±K[<t> + fi\ko)], (2.15) 

where fco is a vector k G A at which /i A T >(k) attains maximum value. 

The joint probability distribution of the random variables a A = {a(r),r£ A} 
is given by the Gibbs measure 

dp^(a A \K,h A ;4) = cxp [-/3H^(a A \K, h A ; *(</>))] [] da(r)/z| T) (^, fe A ; 0), 

feA 

(2-16) 

where da(f) is the Lebesgue measure on R and 

Z A T) (K,h A ;4)= I expU^KI^A;^))] TT da ^ ( 2 ' 17 ) 

is the partition function of the Gaussian model. The latter is finite for all <fi > 
and equals +oo for <fi < 0. The free-energy density of the mean spherical model 
in a finite region A is given by the Legendre transformation 

Pf[ T \K, h A ) := sup {-lAr 1 InZf fc A ;0) - ■ (2.18) 



Here the supremum is attained at the solution <f> — <frT' [K,h A ) (for brevity 

" • ! * "' 



denoted by 4> A ) of the mean spherical constraint 



|A|- 1 ^(a 2 (r*)) A r) (ir,/ lA ^)=l ! (2.19) 



where (• • -)^\k, h\; <f>) denotes expectation value with respect to the measure 

JUS- 

By direct evaluation of the integrals in the partition function (|2.17(l . one 
obtains 

f3f{ T) (K,h A ) = i [ln(K/2n) - K^\k ) 

+U ( A \<P [ A ) ) ~ P ( A \K, h A - - K^}. (2.20) 

Here we have introduced the function 

U-«(^) = |A|- 1 X;M^ + W «(*)] ) (2.21) 
fceA 

which describes the contribution of the spin-spin interaction (called " interaction 
term"), where 

^\k):=ti\k )-ti\k) (2.22) 
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is the normalized excitation spectrum, and the function 

represents the "field term". In (|2.23|) h^\k) denotes the projection of the 
magnetic field configuration ft, a on the eigenfunction (r , k) : 

h£\k) = '£h(r)u ! } [\f,k). (2.24) 

The mean spherical constraint (|2.19|) has the form 

jfPtt) " ^Pi T \K,h A ^) = K (2.25) 

Its solution = <f>T {K, h A ) depends on the lattice region A, the dimensionless 
coupling constant K and the external magnetic field configuration h A . 



3 The finite system with a film geometry 

Hereafter we consider only boundary conditions r = {r\,p, ■ ■ ■ ,p}, with t\ = 
D, N, ND for, respectively, Dirichlct-Dirichlct, Neumann-Neumann and Neumann- 
Diriclct boundary conditions, when ko — {1,^2, •• • ,L c i}. By taking the limit 
^2) • • • » — * oo in expression l|2.21|l at fixed L\ = L we obtain 



lim U (Tl ' p '-- p) ((j)) 



(3.1) 



Next we confine ourselves to the consideration of uniform magnetic fields, h(r) 
h, f € A. By taking the limit of a film geometry in H2.23[l . we obtain 



Pf\K,h;<i>):= T lim P^(K,h A 



1 L 

— Y 



[h^{k)f 



KL t^i <i> + 2 cos <Pl (! ) - 2 cos <^ ( fc ) 



where 



U T \k) := h^u { p{r,k), t £ {D, N, ND}. 

r=l 

From Eqs. (|2.7|) - (I2.11|) we obtain explicitly 



U D \k) = 



L+l 



cot 



irk- 



2(L + 1) 



0, 



/c odd 
k even, 



(3.2) 



(3.3) 



(3.4) 
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and 




TT(fc-l) 

2L 



(-1) 



k-1 



cot 



k = 1 

fc = 2 > .», 
%{2k - 1) 



2(2i + l) 



(3.5) 
(3.6) 



V2T+ 1 

We just mention that for periodic boundary conditions the well-known result 
is P { L p) (K,h;(t)) = h 2 /K(j>. 

Note that due to the field dependence of the solution of the mean spherical 
constraint for the spherical field 0l (K, h) , one has to distinguish between two 
kinds of susceptibilities. The "fluctuation part" of the susceptibility 



t (K,h) 



4=0< T) (i?,/i) 



(3.7) 



measures the fluctuations of the magnetization in the limit of layer geometry 
(more precisely, the variance of a properly normalized block-spin in the corre- 
sponding limit Gibbs state) and satisfies the fluctuation-dissipation theorem. 
On the other hand, by differentiating twice the free energy density with respect 
to the magnetic field, taking into account the implicit dependence on h through 
the solution <f>^[K,K) of the mean spherical constraint, one obtains the total 
magnetic susceptibility per spin, 



k B T X { ll ot (K,h) := -^\pff>(K,h)] = k B T X { lkuJK,h) 



1 d 2 



2d4>dh 



P<?\K,h;4>)] 



dh 



(3.8) 



These two susceptibilities coincide in the zero-field case when there is no spon- 
taneous magnetization in the system. 

For the system with a film geometry the mean spherical constraint l|2.19|l 
takes the form 

W£K4>)-^P£\KM) = K, (3.9) 



where 



1 L 

W £}(<t>) := 7 ^Wd-i[^ + 2cos^)(l)-2cos^ T )(fe)] 



(3.10) 



fc=i 



and 



[•2lT r-2-K 




\ Mi ■ 


■ dd d - x 


'0 Jo 





d-l 



2^(1 - cos e v ) 



v=l 



(3.11) 
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After the evaluation of W^ T J (4>) , the corresponding interaction term t/^ T j (4>^ ) 
in the singular (in the limit L — > oo) part of the free energy density, see (|2.2U|) . 

pf ( :iJ K M = \ [u { :M ] ) - ifHKMP) - k^} , (3.12) 

can be obtained by integration: 

U^M?) = U i: } A<f>°) + l K °W<>)- (3.13) 

Jrj) 

Here <fi^ = <j)^ (K, h) is the solution of equation l|3.9f) , and 4>a > is a suitably 
chosen constant. 

Equations (|3.2|l - (|3.13|l provide the starting expressions for our further finite- 
size scaling analysis. 

4 Critical behavior of a three-dimensional film 

At d = 3 and nonperiodic boundary conditions r = D, N, ND at the surfaces 
of the film, the interaction term (|3.10() in the mean spherical constraint (|3.9|l 
takes the form 

1 L 

W^[<f>) ■= - ^ Wa [<£ + 2 cos ^ T >(1)- 2 cos (4.1) 
fe=i 

This term has been evaluated by using an improved version |l(Jj of the method 
developed by Barber and Fisher [5], |231- Following [5] we set 

W 2 {z) := -(l/47r)lnz+(5/47r)ln2 + Q 2 (z), (4.2) 

where Q2(z) is defined by the above equation. The asymptotic behavior of 
Q2{z) as z — ► follows from the well-known one of the Watson integral W2(z): 

Q 2 (z) = -^-zlnz + 0(z). (4.3) 
Now expression (|4.1|) can be identically rewritten as 

(0) = g[ T) (0) + 9 { 2 T) (0) + (5/47T) In 2, (4.4) 

where 

1 L 

g[ T) W = - 4-r E ln + 2 cos ^i r) (!) - 2 cos ^i T) (*)) ( 4 - 5 ) 
^ k=l 

and 

1 L 

92 ] W = ^ (0 + 2cos^)(l) -2 cos . (4.6) 

fe=l 

We remark that for the boundary conditions under consideration, the func- 
tion (|4.5|l can be calculated exactly, and the function l|4.6l) can be readily evalu- 
ated with the aid of the Poisson summation formula. Here we present the final 
results. 
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4.1 Dirichlet-Dirichlet boundary conditions 

We are interested in the critical regime when cf> — > + and L — ► oo, so that 
0/2 + cos^H 1 ) < L Then we set 



a; = arccos 



-+COS 



L + l 



-,1/2 



(L+l) 



(4.7) 



Under the above substitution, the function g±((p), defined by Eq. 1|4.5|) for 
t = D, reads 



9i 



(D) 



k=l 



2 cos x — 2 cos 



L + l 



(4.8) 



The sum in the right-hand side can be calculated exactly by making use of the 
identity, see [2*5] . 



cos(nx) — cos(ny) = (cos x — cos y) [2 cosx — 2 cos(y + 2fc7r/n)], (4.9) 



fc=i 



setting here ?/ = 0, n = 2(L + 1), and making simple transformations of the 
product with the use of the periodicity of the cosine. Thus one obtains 



1 



4nL 



■In 



sin(L + l)x 



sin x 



(4.10) 



Now we pass to the evaluation of the function g^ D \<f>), defined by Eq. (|4.6[1 
at t = D. Under the substitution (|4.7|l it explicitly reads 



ID' 
92 ' 



1 ( irk \ 

^ = J, z2® 2 ( 2 cos a: - 2 cos J 

2^ 2J ( 2cosx - 2 cos - — ^ j - — Q 2 (2cosx + 2). 



(4.11) 



In deriving the second equality we have used the periodicity of the cosine. For 
4> < n 2 /(L + l) 2 — > + as L — > oo the last term obviously yields Q2(4)/2Z, + 
0(_L -2 ). By applying the Poisson summation formula to the sum in the right- 
hand side of the last equality in Eq. (|4.11|) , changing the integration variable 
and using the periodicity of the integrand, we obtain 



T i -I />7T 

g f\ ( f ) ) = — — / &0Q 2 {2cosx - 2cos( 

L7r J* /(L+l) 



■^[Q 2 (4) - Q 2 (fl] + A 52 D) (0) + 0(L- 2 ), 



(4.12) 
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where 



L + l 



q=1 Jn/(L+X) 



Ltt 



d9cos[2q{L + l)9]Q 2 (2cosx~ 2cos(9). (4.13) 



Consider first the integral in the right-hand side of Eq. (|4.12|) . Its lower 
limit cannot be extended to 0, since we consider the regime when 



2cosz = <j>+ 2cos[7r/(L + 1)] = 2 - x 2 < 2, 



(4.14) 



and the function Q 2 is not defined for negative arguments, see Eq. (|4.2(l . That 
is why we approximate the integrand by 



) 2 (2cosa; — 2cos( 

7T 



Q2 



L + l 



- 2 - 2 cos ( 
2- 2cos6> 



(4.15) 



and notice that the resulting integrals converge when the lower limit tends to 
zero. Then, by expressing Q 2 in terms of W2 through the definition (|4.2() . and 
taking the integral 



1 

4^2 



we obtain 
1 



d9 ln(0 + 2-2cos6>) = In {l + 0/2 +[<£(! + 0/4)] 1/2 } 



1 / 2 

47T 



•0(0). (4.16) 



7T/(L + 1) 



d0Q 2 (2cosa;-2cos0) W 3 (0) 



5 In 2 1 / 2 



4?r 



+ V +0(i } - (417) 



The integral in the right-hand side of Eq. (|4.13|) can be evaluated by twofold 
integration by parts. Since the resulting three terms are oc q~ 2 , the sum over q 
in Eq. (|4.13|l converges and we conclude that 



Agi D) (d>)=0(L- 2 ). 



(4.18) 



By combining Eqs. I|4.12|l . (|4.17|l . I|4.18|l . and taking into account the small 
argument asymptotic form of W3 and Q 2 we obtain 



5 In 2 

47T 



L 



K Ci3 -> 2 (4)-^ 

i <37T 



+ 0(L~ 2 ). (4.19) 



Finally, by substitution of Eqs. I|4.10|l and 14.19|1 into Eq. 14. 411 . and taking 
into account Eq. I|4.7|) , we obtain the zero-field mean spherical constraint in the 
form: 



In 



sin(7r 2 — L 2 ^) 1 / 2 



(tt 2 - L 2 ^) 1 / 2 



-47T 



K, 



c,3 



-W 2 {A) 



InL = 4ttL(K c , 3 - K) 
7 In 2" 



8- 



+ 0(L- 1 ). 



(4.20) 
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This equation coincides (up to a factor of 2, due to the different choice of the 
coupling constant in the Hamiltonian) with the analytical continuation of Eq. 
(7.12) in [5], or Eq. (22) in 9 , from the domain L 2 (j) > ir 2 to the domain 
L 2 (f> < 7r 2 . In that work the focus was on the finite-size scaling behavior and 
the In L term in the left-hand side of Eq. (|4.2U|) was attributed to the finite-size 
shift of the critical coupling: 



InL 1 



AnL L 



(4.21) 



Then, in terms of the variables L 2 (f> and L[K^ L — K) the mean spherical 
constraint (|4.20|) takes the expected finite-size scaling form (in [S] the difference 

K m,l ~ K ' 1S denoted by AK). 

However, Chen and Dohm 12^ have noticed that at the bulk critical temper- 
ature, when K = K c3l Eq. (|4.20J) has a leading-order solution <f>^ {K c ^, 0) ~ 
L~ 3 . This behavior follows by assuming L 2 (f> — > as L ~ * oo and expanding the 
left-hand side of the mean spherical constraint up to the leading order: 

\n<j) + 31nL = 4nL(K c , 3 - K) + 0(1). (4.22) 

Note that, if K c3 > K, i.e. L(K C ^ — K) — » oo, the solution of the above 
equation is 

0~ L -3 e -4irL(tf c , 3 -x)+o(i). (4.23) 

To relate the above fact to the critical behavior of the model, we start from 
the exact result for the zero-field susceptibility per spin under Dirichlct-Dirichlct 
boundary conditions, see Eq. IjA.lOf) . 



y Wt K 0) = — r cot 2 W2)sin[(£ + l)x] _ L + l , 
Xl [ ' > 2j\Lsmx[l + cos(L + l)x} 2L s\n 2 {x/2) ' ' '"" 



By expanding the right-hand side for L cf> — ► as L — ► oo, we obtain 

X ( L D \K, 0)^-^ + ((^r 1 ) + O (L- 2 ) . (4.25) 
Hence, at the bulk critical temperature one recovers the leading-order result of 

B21 

X { L D) (K c , 3 ,0)~L\ (4.26) 

in a full accordance with the predictions of the phenomenological finite-size scal- 
ing, see Eq. ijTTPo)) . For K Ct3 > K one has (K < K c , 3 , 0) - L 3 exp[47ri(^ Ci3 - 
K)], i.e. the susceptibility diverges then exponentially, which agrees with Eq. 
(| 1 . 1 311 in the case d! = 2. 

To understand this anomalous behavior at T c , we recall that the "standard" 
finite-size scaling regime takes place when L(K^' L — K) = 0(1), or equivalently, 
at 

K = K C , 3 - 1 ^-^, Xl =0(l), (4.27) 
4ttL L 
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i.e., in a narrow temperature interval of width 0(L _1 ) (note that in the three- 
dimensional spherical model v = l/(d — 2) = 1) which is shifted by 0{\aL/L) 
above the bulk critical temperature T c3 . Thus, the change to the behavior 
Xl ~ L 3 , which takes place at the bulk critical point, occurs in a temperature 
region of the same width shifted by 0(lnL/L) below the finite-size scaling one. 

We emphasize that the result Xl (K c3 , 0) ~ I? is characteristic of the layer 
geometry with two free surfaces. To show that we consider below the cases of 
Neumann-Dirichlet and Neumann-Neumann boundary conditions. 



4.2 Neumann-Dirichlet boundary conditions 



In the critical regime when < 


f> -> 0+ and L - 


■+ CO, 


we set 












( n V 




x — arccos 




- COS 








\2L+lJ_ 





(2L + 1) 2 



1/2 



(4.28) 



In this case the exact expression for the function g[ ND \(t>), defined by Eq. 
H4.5JI . follows from the identity (4.20) in [TO] : 



4nL 



■In 



cos(L + l/2)x 
cos(x/2) 



(4.29) 



where x is given by Eq. I|4.28|l . The evaluation of g^ D ^ (<fi) goes along the same 
lines as in the previous case with the result 



1^(0) = *ca 



3 In 2 1 
~4^T + 2L 



K, 



c,3 



In 2 
"2V 



0{L~ 



(4.30) 



Now the zero-field mean spherical constraint takes the form (see Eq. (4.22) in 
[TO] at zero surface fields and replace there 2K — > K) 



In 



cos(tt 2 /4 - L 2 <j)) 1/2 \ = 4ttL(K c>3 - K) + 2tt 



K, 



c,3 



W 2 (A) 



ln2" 

OiL-HnL). (4.31) 



Comparing the above equation with the analogous Eq. (|4.2U|) . we see that 
now there is no hiL/L finite-size shift of the critical temperature: 



K 



(ND) _ 



K c , 3 -W2(4)- 1 ^ 

ZTT 



(4.32) 



The mean spherical constraint i|4.31ll takes the expected finite-size scaling form 
in terms of the variables L 2 (j) and L{K^^ — K). Even in the regime L 2 <p — > 
as L — > 00, by expanding its left-hand side one obtains in the leading order 



ln< 



21nL = AttL{K c ^ - K) + 0(1). 



(4.33) 
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Therefore, at the bulk critical temperature one has the standard finite-size be- 
havior (if c ,3) 0) ~ L~ 2 . Note that below T c , when L{K C ^ — K) — > oo, the 
solution of the spherical field equations is 



L -2 e -4irL(K C:3 -K)+0(l) 



(4.34) 



i.e. it is similar (only the power of L in front of the exponential term differs) to 
that of Dirichlet-Dirichlet boundary conditions, see Eq. I|4.23[) . 

The critical behavior of the zero-field susceptibility per spin under Neumann- 
Dirichlet boundary conditions follows from the exact result, see Eq. (IA.19|I . 



4cot 2 (x/2) 



2H2L+ l)cos 2 (a;/2) 



tan[(L + l/2)a;] 
Lsmx 
L + l 



2(2L + l)sin^(a;/2) 



(4.35) 



where x is given by Eq. (|4.28(l . By expanding the right-hand side of the above 
equation for a; < 1, we obtain to the leading order 



J(tt 2 /4-L 2 



tan( 7 r 2 /4-L 2 0) 1 / 2 
(tt 2 /4- L 2 ^) 1 / 2 



1 



(4.36) 



From Eq. I|4.31(l it is clear that for any finite value of L(K C ^ — K) the suscepti- 
bility diverges as L 2 when L — >• oo. Note that (tt 2 /4 - L 2 ^) 1 / 2 -> 0+ only when 
L(K^l — K) — ► 0. Therefore, at the shifted critical temperature one obtains 
the simple result 

L 2 



(ND),Tf 



(ND) 
rn , L ' 



o)- 



3J 



(4.37) 



Thus we conclude that in the presence of only one free (Dirichlet) surface the 
mean spherical model has the usual finite-size critical behavior. It is instruc- 
tive to see whether this behavior will change in the presence of two equivalent 
surfaces with Neumann boundary conditions. 



4.3 Neumann-Neumann boundary conditions 

The interaction term l|4.1|l in this case has been treated completely and rigor- 
ously in [Till. Since cos 1/3^(1) = 1, see (|2.11|) . we set 



cosh -1 (l 



a/a 



(4.38) 



where cosh 1 denotes the inverse function of cosh. The exact expression for the 
function g[ N \(j)) is given by Eq. (3.13) of ^U]. It reads 



9i (4>) 



4ttL 



ln< 



In 



sinh Lx 
sinhx 



(4.39) 
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In the critical regime <j> — ► as L — > oo the result for the function g% (0) i s 
(see Eq. (3.18) in JOj under the replacement 2K C — ► if c ): 



if. 



5 In 2 

47T 



1 

2L 



W a (4) 



3 In 2 

47T 



0(L" 



(4.40) 



Thus, by substitution of equations 1)4.39)1 and (|4.4Q|I into 1)4.4(1 . in the limit 
— * 0, L — > oo, we obtain in the leading order the following finite-size scaling 
form of the zero- field mean spherical constraint: 



In 



L^l 2 sinh(I^ 1 ' a ) 



1/2^ 



where K^' L is the shifted critical coupling, 



K 



(AT) 



lni 1 
Kc ' 3 + ^L~2L 



Wa(4) 



3 In 2 



(4.41) 



(4.42) 



Comparing Eqs. 1)4.42)1 and 1)4.2111 we see that under Neumann-Neumann 
boundary conditions the finite-size scaling region of width 0(L _1 ) is shifted 
below the bulk critical temperature T Cj 3 by Oi\nL / L). Since at the bulk critical 
point, when K = -Kc.3, the right-hand side of Eq. 1)4.41)1 goes to plus infinity 
like InL as L — > oo, we conclude that L^ 1 / 2 — > oo, although <f> — ► + . In this 
regime Eq. 14.4111 simplifies to 



]n^ 2 + L^ 2 = 0(l), 



which yields 



{\nL/Lf 



(4.43) 



(4.44) 



This, at the bulk critical temperature, the finite-size critical behavior of the 
spherical field under Neumann-Neumann boundary conditions becomes loga- 
rithmically modified. Note that below T c , i.e. when L(K c> a — K) — > oo, the 
solution of the spherical field equation (|4.41)l is 



L -l e -4nL(K , 3 -K)+0(l) 



(4.45) 



i.e. the leading-order behavior of <f> again differs only in the power of L in front 
of the exponential term, compare with Eqs. 1)4.23)1 and 1)4.34)1 . 

The corresponding critical behavior of the zero-field susceptibility follows 
from the exact expression, see (|A.27f) . 



x£°(*,o) 



1 1 + cosh : 



L 2 sinh : 



[L coth(Lir) — cotha;] 



(4.46) 



where x is given by Eq. 1)4.38)1 . By expanding this expression for x = (j) 1 / 2 -> 0+ 
and Lx = Lcj) 1 ^ 2 — ► oo, we obtain the following modified leading-order behavior 
of the finite-size zero-field susceptibility at the bulk critical point: 



xi?°(*c.,0)S 



1 



L 



J(lnL) 



(4.47) 
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in a full accordance with Eq. (|1.12|) . 

We emphasize that the above behavior takes place at the bulk critical tem- 
perature, in a temperature region of width 0(L _1 ) shifted by (9(lnL/L) above 
the finite-size scaling one. It will be further modified, although not radically, 
if one considers the temperature interval located at p > times the same shift 
above T C) 3. Then one obtains from Eq. Q4.41fl . up to the leading order, 

<p { p{K L , p )~[(\+p)\nL/L]\ (4.48) 

and hence the corresponding reduction of the finite-size zero-field susceptibility 
(|4"47|) by a factor of (1 +p)~ 2 . 

5 Discussion 

In the current article we have presented the predictions of the phenomenological 
finite-size scaling for systems with asymptotically large shift of the bulk critical 
temperature. We have shown that in such systems the behavior of the zero- 
field susceptibility at the bulk critical point depends crucially on the sign of the 
shift - the positive shift leads to the reduction of the "standard" divergence of 
~ L 1 ^ , while the negative shift leads to a stronger divergence. We have 
verified our considerations on the example of the three-dimensional spherical 
model under Dirichlct-Dirichlct, Dirichlct-Ncumann, and Neumann-Neumann 
boundary conditions. A Dirichlet surface leads to a negative shift of the critical 
coupling, — lnL/(87rL), and a Neumann surface - to a positive one, lnL/(87rL) 
jllj . So, for a film geometry under Dirichlet-Dirichlet boundary condition the 

shift 4 D) = -]nL/(4nL) leads to x ( l\k c ^Q) ~ L 3 , (see Eq. KZty . in full 
accordance with our phenomenological predictions (see Eqs. I|1.13(l and Ijl.lSp V 
For Neumann-Neumann boundary conditions the shift is positive, i.e. = 

lnL/(47r_L), wherefrom Xl (T c ) — L 2 / '[J (In L) 2 ] , in a full agreement with the 
phenomenoligal prediction given by Eq. (|1.12(l . 

We emphasize that, in order to make concrete phenomenological prediction 
for the finite-size behavior of the zero-field susceptibility in systems with asymp- 
totically large shift of the critical temperature, use has been made of the the 
size dependence of Xl f° r T <T C . We have supposed that the leading behavior 
of x^ T ' will be the same under other boundary conditions. We have verified this 
assumption on the example of the spherical model. It is highly desirable to have 
the corresponding results for other models too. 

One might ask why is the shift of the critical coupling in the spherical 
model positive under Neumann-Neumann and negative under Dirichlct-Dirichlct 
boundary conditions. Indeed, this contradicts the general expectations based 
on arguments like that the missing neighbors in a ferromagnetic system should 
reduce its critical temperature. In order to understand the above facts, let us 
note that the observed behavior is in agreement with the length of the spins near 
the boundary. This length is reduced near a Dirichlet boundary and enlarged 
near a Neumann one (there < a 2 >~ 1.34 [23]; similar estimation for the 
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Dirichlet boundary gives < a 2 >~ 0.83). Then, since the total length of all the 
spins is fixed, that leads to spins in the main part of the system being larger 
than 1 under Dirichlet and smaller than 1 under Neumann boundary conditions. 
As a result, an effective interaction is taking place with spins which length is 
not equal to 1. In turn, this produces a shift in the "critical temperature" of the 
finite system which is positive for Dirichlet and negative for Neumann boundary 
conditions, contrary to what one would expect for a system with a fixed length 
of the spins. 

Finally, let us recall that the infinite translational invariant spherical model 
is equivalent to the n — > oo limit of the corresponding system of n-component 
vectors [2|)],[2Z] 5 but the spherical model with surfaces (or, more generally, with- 
out translation-invariant symmetry) is in fact not such a limit |28| . In other 
words, the spherical model under nonperiodic boundary conditions is not in the 
same surface universality class as the corresponding O(n) model in the limit 
n — » oo, in contrast with the bulk universality classes. The last becomes appar- 
ent when one investigates surface phase transitions for an 0(n) model in the 
limit n — > oo. In that case in the limit n — > oo one obtains (55] A° = l/(d — 2) 
(i.e. A? = 1 for d = 3) for ordinary and Af b = 2/(d - 2) (i.e. Af = 2 for 
d = 3) for special phase transitions, while A° = 1/2 and A| b = 3/2 ^U], [231 
for the three-dimensional spherical model. It is believed that the corresponding 
equivalence will be recovered if one imposes spherical constraints in a way which 
ensures that the mean square value of each spin of the system is the same |2*5) 
(unfortunately such a model is rather untractable analytically). The introduc- 
tion of just one additional spherical field that fixes the mean square values of 
the spins at the boundaries changes the surface critical exponents but is 
not enough to recover the correspondence to the 0(n) models. Unfortunately, 
the finite size scaling properties even of that analytically tractable model have 
not been investigated. 
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A Calculation of the field term 

A.l Dirichlet-Dirichlet boundary conditions 

In uniform magnetic field h and in the regime when cos a; := 0/2 + cosy^ (1) < 
1, the field term Eq. 13.211 has the explicit form (assuming L odd) 



, 2 (L-l)/2 .2 7r(2fc+l) 
COt 2 (L + 1) 



KL(L+1) 



IT £ 



k=0 



COS x — COS 



r(2fc+l) 
L+l 



By making use of the elementary identity 

cot 2 (a/2) = 



1 — cos a 

the above expression can be rewritten as 



1, 



(A.l) 



(A.2) 



where 



K 1 1 — cos x 



S?\x)-S?\0)]-S?\ X ) , (A.3) 



(L-l)/2 



— r 

v + 1 A: 



i ( L + 1 ) fe=0 cos x cos ■ j^rp 



7r(2fc+l) 



(A.4) 



To calculate the sum (x) , we use the identity (|4.9|) and set there n = L + 1 



and y = n/(L + 1). This yields 



n 

k=0 



cos x — COS ■ 



7r(2fc + 1) 



L + l 



2" i [cos( J L + l)x + l]. 



(A.5) 



Since the left-hand side of the above identity equals 



(L-l)/2 

n 

fc=0 



COS x — COS 



7r(2fc + 1) 



L + l 



(A.6) 



we obtain 



(L-l)/2 

n 

fe=0 



cos x — cos ■ 



7r(2fc + 1) 



L + l 



2- L/2 [cos(L + l)x + l] 1/2 . 



(A.7) 



21 



Next, by taking logarithm of both sides of Eq. (|A.7|I and differentiating the 
result with respect to x, we obtain 



sin(L + l)x 



Hence 



S^(0) = BmSr(*) = =^ 



2Lsina;[l + cos(£ + l)x] ' 



(A.. 



(A.9) 



Finally, by inserting Eqs. (| A. 8|) and ljA.91) into Eq. (|A.3|) . and making elemen- 
tary transformations based on the identity (| A. 2|) . we obtain the exact result 



cot 2 (a;/2)sin[(L + l)x] 



L + l 



(A.10) 



2K \Lsina:[l + cos(L + l)a;] 2isin 2 (a;/2) 

A. 2 Neumann-Dirichlet boundary conditions 

In uniform magnetic field h and in the regime when cos x := 0/2+cos tp^ (1) < 
1, the field term Eq. I|3.2[l has the following explicit form 



P[ ND \K,h;4>) = 



L-l 

E 



cot 



2 7r(2fc+l) 
2(2L+1) 



^(2i+l)^ cos ,_ cos ^ 



(A.ll) 



By making use of the elementary identity (| A. 2|> the above expression can be 
rewritten as 

h 2 



p( ND \K,h-< 
where 

5' 



K I 1 — cos x 



S 



(ND) 



{x)-S[ ND \0) 



s 



(ND) 



(X) 



iND \x) = 



1 

IL + I) E ' 



7T(2fc+l) 



(A.12) 
(A.13) 



L w 2L(2L+ 1) ^ cos a; -cos- 2L+1 

To calculate the sum S^ D \x), we use the identity (|4.9|) and set there n = 2L+1 
and y = ir/(2L + 1). This yields 

2L 



n 



fc=0 



COS x — COS 



7r(2fc + 1) 



2" ijL [cos(2L+ l)x+ 1]. 



2L + 1 

Since the left-hand side of the above identity equals 

(L-l 



(cosx + 1) <^ J| 



. fe=0 



cos x — cos ■ 



7r(2fc+ 1) 



2L+ 1 



we obtain 

L-l 



n 



cos a; — cos • 



7r(2fc + 1) 



2L+ 1 



>-L 



cos(2L + l)x + 1 



cos x + 1 



1/2 



(A.14) 



(A.15) 



(A.16) 
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Next, by taking logarithm of both sides of Eq. I|A.16|) and differentiating the 
result with respect to x, we obtain 



S 



Hence 



(ND) 



(,:) 



sin(2L + l)x 



1 



4Lsina;[l + cos(2L + l)x] AL(2L + 1)(1 + cos: 



S 



(ND) 



(0) = ]imS ( T ND \x) = n , 
w x->o L w 2(2L + 1) 



L + l 



(A.17) 



(A.18) 



Finally, by inserting Eqs. (|A.17|) and i|A.19|l into Eq. (|A. 12|) . and making 
elementary transformations, we obtain the exact result 



pW{K,h>, 



K 



4cot 2 (x/2) 



2L(2L + l)cos 2 (x/2) 



tan[(£ + l/2)af] 
Lsinx 
L + l 



2{2L+ l)sirr(x/2) 



(A.19) 



A. 3 Neumann-Neumann boundary conditions 

This case differs from the previous two ones in that now cos (1) = 1 and we 
have to set cosh a; := 1 + 4>/2 > 1. In uniform magnetic field h, the field term 
Eq. 11 -J. 21 has the explicit form 



h 2 



2h 2 



K<j> K 



where 



fc=i 



cos 



2 TT_k 
2L 



coshx — cos(7rfc/L) 



By making use of the elementary identity 

cos 2 (a/2) = -(1 + cosh x) — — (cosh x — cos a) 
the above expression can be rewritten as 

L-l 



(A0 _ 1 + COshx ^ 
S L W ~ 2L 2 



1 



fc=l 



L — 1 



cosh a; — cos(nk/L) 2L 2 



To calculate the above sum, we start from the identity 25 



(A.20) 



(A.21) 



(A.22) 



(A.23) 



n 



fc=i 



z 2 — 2z cos 



77 k 



y 2n 



z 2 -l 



(A.24) 



By setting here z — exp(x), n — L, and performing elementary transformations, 
we obtain 

coshx-cos— =2~ L+1 ^-- — '-. A.25 
L sinn x 



n 

fc=i 



23 



Next, by taking logarithm of both sides of Eq. (|A.7|I and differentiating the 
result with respect to x, we obtain 

r,(N) f \ 1 + cosh x . , . , , L — 1 

5 ' i = — — — \L coth (is) - cothx - — — . (A.26 
L v ; 2L 2 smha; L v ' 1 2L 2 y J 

Finally, by inserting equation l|A.26|l into equation 1|A.20|) . we obtain the exact 
result 

Pf\K M ) = ^{\ + [icoth(Lx) - coth,] - ^1} . (A.27) 



24 



